Fidelity approach to quantum phase transitions: finite size scaling for quantum Ising model in a 
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We analyze the scaling parameter, extracted from the fidelity for two different ground states, for the one- 
dimensional quantum Ising model in a transverse field near the critical point. It is found that, in the thermo- 
dynamic limit, the scaling parameter is singular, and the derivative of its logarithmic function with respect to 
the transverse field strength is logarithmically divergent at the critical point. The scaling behavior is confirmed 
numerically by performing a finite size scaling analysis for systems of different sizes, consistent with the con- 
formal invariance at the critical point. This allows us to extract the correlation length critical exponent, which 
turns out to be universal in the sense that the correlation length critical exponent does not depend on either the 
anisotropic parameter or the transverse field strength. 
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Introduction. An emerging picture arises due to latest ad- 
vances in quantum information science, which allows us to 
study quantum phase transitions (QPTs) [ 1 ] from the ground 
state wave functions of many-body systems. One of the well- 
studied aspects is to unveil the possible role of entanglement 
in characterizing QPTs H H H fy, H @] (for a review, see fl). 
Remarkably, for quantum spin chains, the von Neumann en- 
tropy, as a bipartite entanglement measure, exhibits qualita- 
tively different behaviors at and off criticality J3l . 

On the other hand, the fidelity, another basic notion of 
quantum information science, has attracted a lot of atten- 
tion fgl [Tol EH quite recently. In Ref. ifioll . it has been shown 
that it may be used to characterize QPTs, which occur in 
quantum spin chain, regardless of what type of internal order 
is present in quantum many-body states (either the conven- 
tional symmetry-broken orders or exotic QPTs in matrix prod- 
uct systems [12]). The argument is solely based on the basic 
Postulate of Quantum Mechanics on quantum measurements. 
Indeed, the basic Postulate of Quantum Mechanics on quan- 
tum measurements implies that two non-orthogonal quantum 
states are not reliably distinguishable II 1311 . Therefore, any two 
ground states must be orthogonal due to the occurrence of or- 
ders, regardless of what type of QPTs. Conversely, the fact 
that two ground states are orthogonal implies that they are re- 
liably distinguishable. Therefore, an order parameter, which 
may be constructed systematically in principle 114I1 . exists for 
any systems undergoing QPTs. It is the quantitative or quali- 
tative difference unveiled in order parameters that justifies the 
introduction of the notions of irrelevant and relevant informa- 
tion. To quantify irrelevant and relevant information, the scal- 
ing parameter extracted from the fidelity was introduced to 
characterize QPTs. This establishes an intriguing connection 
between quantum information theory, QPTs, renormalization 
group (RG) flows and condensed matter physics. 

The fact that any two different ground states are orthogonal 
for continuous QPTs makes it difficult (if not impossible) to 
extract physical information solely from ground states them- 
selves. Conventionally, condensed matter physicists and field 



theorists focus on spectra and correlation functions. There- 
fore, it is somewhat surprising to see that simply partitioning 
a system into two parts and quantifying entanglement between 
them reveal highly nontrivial information about QPTs. The in- 
trinsic irreversibility due to information loss along RG flows 
may also be revealed solely from ground states I4l Il5lll6l.il 711 . 



In the fidelity approach 11011 . it is necessary to put the whole 
system on a finite chain, and observe how the fidelity scales 
with system sizes as the thermodynamic limit is approached, 
in order to extract physical information. The difference be- 
tween entanglement measures and the fidelity approach lies in 
the fact that for the former different entanglement measures 
need to be devised to detect QPTs H, whereas the latter suc- 
ceeds to detect QPTs for quantum spin chains, regardless of 
what order is present. The philosophy behind this is that bi- 
partite entanglement measures involve partitions and some in- 
formation is lost due to the fact that the whole is not simply the 
sum of the parts, whereas in the fidelity approach, a system is 
treated as a whole from the starting point. 



In this paper, we analyze the scaling parameter, extracted 
from the fidelity, for the one-dimensional quantum Ising 
model in a transverse field near the critical point. It is found 
that, in the thermodynamic limit, the scaling parameter is sin- 
gular, and the derivative of its logarithmic function with re- 
spect to the transverse field strength (the control parameter) 
is logarithmically divergent at the critical point. A finite size 
scaling analysis is carried out for systems of different sizes, 
and the scaling behavior is confirmed numerically, consistent 
with the conformal invariance at the critical point. This allows 
us to extract the correlation length critical exponent. We have 
also performed numerics to confirm the universality, i.e., the 
correlation length critical exponent does not depend on either 
the anisotropic parameter or the transverse field strength. 



The fidelity and the scaling parameter for quantum XY spin 
chain. The quantum XY spin chain is described by the Hamil- 



2 



tonian 



H ■■ 



M 

j=-M 



1+7 



(1) 



Here cr x ,cr., and cr 2 . are the Pauli matrices at the /'-th lat- 
tice site. The parameter y denotes an anisotropy in the 
nearest-neighbor spin-spin interaction, whereas A is an ex- 
ternal magnetic field. The Hamiltonian (HJ may be exactly 
diagonalized JH, 19] as H — % k A k (c k c k - 1), where = 
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(A - cos(27rk/L)) 2 + y 2 sin ilnkjL), with c k and c\ denot- 
ing free fermionic modes and L = 2M + 1 . The ground state 
is the vacuum of all fermionic modes defined by c k \if/) = 
0, and may be written as \if/) = T\¥=i(cos(0 k /2)\0) k \0)- k - 
i s,in(6 k /2)\l) k \l)- k , where \Q) k and \l) k are, respectively, the 
vacuum and single excitation of the k-th mode, and 8 k is de- 
fined by cos k = (cos(2j:k / L) - A)/A k . Therefore, the fidelity 
F for two different ground states \if/(A, y)) and \if/(A', y)) takes 
the form: 



u 



F{A,X;y) = \\ 



cos ■ 



(2) 



where the prime denotes that the corresponding variables take 
their values at A'. Obviously, F = 1 if A = X. Generically, 
cos -^s-i < 1, therefore the fidelity (O decays very fast when 
A separates from A'. 

Now let us introduce a fundamental quantity-the scaling pa- 
rameter d(A, A'; y). For a large but finite L, the fidelity scales 
as d L , with some scaling parameter d depending on A and A', 
due to the symmetry under translation. Formally, in the ther- 
modynamic limit, d(A, A') may be defined as 



\nd{A,A'\y) = lim \nF(A,X;y)/L. 



(3) 



The scaling parameter d(A, A'; y) enjoys some properties in- 
herited from the fidelity: (1) symmetry under interchange 
A < — > X; (2) d(A,A;y) = l;and (3) < d(A, A'; y) < 1. 

In the thermodynamic limit, the scaling parameter 
d(A, A'; y) for the quantum XY model takes the form: 

1 r 

\nd(A,A';y) = — dalnT(A,A';y;a), (4) 
2tt Jo 



where 



T(A, A'; y; a) = cos[§(A; y; a) - §(A'\ y; a)]/2, 



(5) 



with 



cos &(A; y; a) = (cos a - A) I -\j(cos a - A) 2 + y 2 sin 2 a. (6) 

A notable feature of the scaling parameter (0]l is that, besides 
d(A,A';y) = d(A',A;y) and d{A,A;y) = 1, it even detects the 
duality between two phases A > 1 and A < 1 for quantum 
Ising model in a transverse field (y — 1) [ 19], since it satisfies 
d(A,A';l) =d(l/A,l/X;l). 



It has been shown ifioll that the scaling parameter d(A, A'; y) 
exhibits a pinch point at (1, 1), i.e., an intersection of two sin- 
gular lines A — 1 and X = 1, for quantum Ising model in a 
transverse field (y = 1). In Fig. [2 we plot the scaling param- 
eter d(A, A'\y) against A for different values of X and y. One 
observes the continuity, as it should be for continuous QPTs. 




FIG. 1: (color online) The scaling parameter d(A,A';y), extracted 
from the fidelity for two ground states \4i(X)) and \ifr{A')) of quantum 
Ising model in a transverse field, is regarded as a function of A for 
some fixed values of A' and y. It is continuous but not analytic at 
A c = 1. (a): the red line is for A' = 2,y = 1, which touches the 
blue line at A = 2 and the green line is for X = 1/2, y = 1, touching 
the blue line at A = 1 /2. The mirror symmetry between two curves 
results from the duality, (b): the green line is for A' = 1/2, y = 1 /2, 
touching the blue line at A = 1/2 and the red line is for X = 2, y = 
1/2, touching the blue line at A = 2. No mirror symmetry for y + 1. 

Let us focus on the quantum Ising universality class with 
the critical line y + and A c = 1. There is only one (second- 
order) critical point A c = 1 separating two gapful phases: (spin 
reversal) Z2 symmetry-breaking and symmetric phases. The 
order parameter, i.e., magnetization (cr x ) is non-zero for A < 
1, and otherwise zero. At the critical point, the correlation 
length ^ ~ \A - A c \ v with v — 1 [19:]- Our purpose is to extract 
the correlation length critical exponent by performing a finite 
size scaling analysis for d(A, A'; y). 

Finite size scaling. In order to quantify the drastic change 
of the ground state wave functions when the system undergoes 
a QPT at the critical point A c = 1 , we evaluate the derivative 
of In d{A, X;y) with respect to A. In the thermodynamic limit, 
]nd(A, X;y) is logarithmically divergent at the critical point 
A r = 1: 



d\nd(A, X;y) 
OA 



k\ In \A - A c \ + constant, 



(7) 



where the prefactor k\ is non-universal in the sense that it 
depends on X and y. The numerical results are plotted in 
Fig- ID for X = 2 and y = 1. The least square method yields 
k\ x, -0.079742. For systems of finite sizes L's, there are no 
divergence in the derivatives of In d(A,A';y) with respect to 
A, since the second-order QPT only occurs in the thermody- 
namic limit. Instead, some pronounced peaks occur at the so- 
called quasi-critical points A m that approach the critical value 
as A m ~ 1 — 5.52233L -0 " 321 , with the peak values logarith- 
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FIG. 2: (color online) Main: the logarithmic divergence near the 
critical point X c = 1 is analyzed. This is achieved by consider- 
ing d A lnd(A,A' = 2;y = 1) as the function of the transverse field 
strength A. The curves shown correspond to different lattice sizes 
L = 201,401,1201,2001,4001,oo. The maximum gets more pro- 
nounced, with the system size increasing. Inset: the position of 
maximum approaches the critical point Ac = 1 as A m ~ 1 - 
5.52233L- " 321 . 

mically diverging with increasing system size L, 



d\nd(A, A';y) 



OA 



\a=a, 



— £2 In L + constant, 



(8) 



where the non-universal prefactor ki takes the value ki ~ 
0.079773. The scaling ansatz in the systems exhibiting log- 
arithmic divergences [20] requires that the absolute value of 
the ratio k\lk 2 is the correlation length critical exponent v. In 
this case, \k\lk2\ ~ 0.999613, very close to the exact value 1. 

In the case of logarithmic divergences, a proper scaling 
ansatz has been addressed in Ref. 12011 . Taking into ac- 
count the distance of the maximum of d^\nd{A,A';y) from 
the critical point, we choose to plot 1 - exp[(9,i In d(A, A'; y) - 
d^hidiA, /t';y)U=/i m ] as a function of L(A - A m ) for different 
system sizes L's. All the data for different L's collapse onto a 
single curve. The numerical results for the size ranging form 
L = 201 up to L = 4001 are plotted in Fig. [3] All these indi- 
cates that the system is scaling invariant, i.e., |/L = (and 
thus conformally invariant), and that the correlation length 
critical exponent v = 1 . 

Universality. As is well known, the quantum XY chain be- 
longs to the same quantum Ising universality class for non- 
zero y, with the same critical exponents. To confirm the uni- 
versality, we need to check the scaling behaviors for differ- 
ent values of y. For A' — 2 and y = 1/2, in the thermo- 
dynamic limit, it takes the form (0 with k\ * -0.157162, 
as long as the control parameter is close to the critical point, 
whereas for a system of finite size, it takes the form (0 with 
&2 ~ 0.157176. Thus, the absolute value of the ratio fci/fo is 
\ki/k 2 \ = 0.999910. Fig.fflshows that all the data for different 
L's collapse onto a single curve. We also plot the derivative 
of the logarithmic function of the scaling parameter d(A, A';y) 
with respect to A for A' — 2 and y = 1/2 (see the inset in 



FIG. 3: (color online) Main: a finite size scaling analysis is carried 
out for a quantity defined as D(A, A';y) = 1 - exp[d A In d(A, A'\y) - 
d A lnd(A,A'\y)\ A=Aw ]. According to the finite size scaling ansatz in 
the case of logarithmic divergences, one expects that D(A,A'\y) is 
a function of L(A - A m ). Indeed, all the data from L = 801 up to 
L = 4001 collapse on a single curve. This shows that the system at 
the critical point is scaling invariant (and thus comformally invariant) 
and that the correlation length critical exponent v is 1. Inset: the 
peak value of d A \nd{A,A' = 2;y = 1) at A,„ diverges as the system 
size increases, leading to k 2 ~ 0.0797726. 
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FIG. 4: (color online) The universality hypothesis for the scaling 
parameter extracted from for the quantum Ising model in a transverse 
field is checked against different values of y and A'. Main: in this case 
we have chosen y = 1/2 and L ranging from 2801 up to 6001. All 
the data collapse, consistent with the fact that the correlation length 
critical exponent v is 1. The inset shows that the derivative of the 
logarithmic function of the scaling parameter d(A, A';y) with respect 
to A for A' = 2 and y = 1/2 is logarithmically divergent at A c = 1, 
with A m ~ 1 - 3.23906ZT 1 ° 1135 . 



Fig.|4|. All the above results show that the critical exponent 
v= 1. 

Besides y, the scaling parameter d(A, A'; y) also depends on 
the control parameter A'. For A' = 1/2 and y = 1, in the ther- 
modynamic limit, the derivative of the logarithmic function of 
the scaling parameter d(A,A';y) with respect to A still takes 
the form (f7j), with k\ w 0.083005, as long as the control pa- 
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rameter is close to the critical point, whereas for a system of 
finite size, it takes the form ©, with k 2 ~ -0.083007. Thus, 
the absolute value of the ratio k 2 /h is \k 2 /h\ = 0.999975 , 
again close to the exact value 1 . Similarly, all the data for dif- 
ferent L's collapse onto a single curve, as shown in Fig. [5] In 
the inset, we plot the derivative of the logarithmic function of 
the scaling parameter d(A, A'; y) with respect to A for A' - 1/2 
and 7=1. Therefore, we have demonstrated that the univer- 
sality hypothesis is valid for the scaling parameter d(A, A'; y). 



o.o 




for the quantum Ising model in a transverse field. In this re- 
gard, algorithms for periodic boundary conditions [22] and in- 
finite systems 12311 are powerful enough to extract meaningful 
information for critical systems. 

In summary, we have performed a finite size scaling analy- 
sis for the scaling parameter, whose analytical expression has 
been extracted from the fidelity for two ground states corre- 
sponding to different values of the control parameter for the 
one-dimensional quantum Ising model in a transverse field 
near the critical point. In the thermodynamic limit, the loga- 
rithmical divergence of the derivative of the scaling parameter 
with respect to the transverse field strength is demonstrated 
numerically, consistent with the conformal invariance at the 
critical point. This makes it possible to extract the correlation 
length critical exponent. The latter turns out to be universal, 
in the sense that the correlation length critical exponent thus 
extracted does not depend on either the anisotropic parameter 
y or the transverse field strength A. 

We thank Sam Young Cho, John Fjaerestad and Jon Links 
for helpful discussions and comments. 



FIG. 5: (color online) The universality for the scaling parameter is 
checked against different values of y and A' . Main: in this case we 
have chosen y = 1 and L ranging from 2801 up to 6001. Consistent 
with the universality hypothesis for the quantum Ising model in a 
transverse field, all the data collapse, indicating that the correlation 
length critical exponent v is 1. The inset shows that the derivative 
of the logarithmic function of the scaling parameter d(A,A';y) with 
respect to A for A' = 1/2 and y = 1 is logarithmically divergent at 
A c = 1, with,!,,, ~ 1 + 3.50186ZT 94107 . 

Discussions and conclusions. As a basic notion of quan- 
tum information science, fidelity may be used to detect QPTs 
in condensed matter systems. Remarkably, an intimate con- 
nection exists between RG flows, QPTs and the scaling pa- 
rameter which may be extracted from the fidelity II X Ofl . The 
scaling parameter is well defined in the thermodynamic limit, 
in sharp contrast to the fidelity itself that always vanishes for 
continuous QPTs. Different from a bipartite entanglement 
measure, the fidelity approach does not involve the partition of 
the whole system into different parts, and the system is treated 
as a whole from the starting point. In some sense, such a dif- 
ference may be counted as the contribution from multipartite 
entanglement. Therefore one may expect that the fidelity ap- 
proach possesses significant advantage over the conventional 
bipartite entanglement approach 12111 . 

Another feature worth to be mentioned is that fidelity is 
simple to be evaluated in the matrix product state (MPS) rep- 
resentation [10]. On the other hand, many efficient numerical 
algorithms are now available due to the latest developments 
in classical simulation of quantum systems 1E2I I23I 12411 . This 
makes it practical to determine all information including sta- 
ble and unstable fixed points along RG flows [10], and to ex- 
tract critical exponents from the scaling parameter, as shown 
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